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Abstract 

Quantum gravity coupled to scalar massive matter fields is investigated 
in the framework of causal perturbation theory. One-loop calculations in- 
clude matter loop graviton self-energy and matter self-energy and yield 
ultraviolet finite and cutoff-free expressions. Perturbative gauge invariance 
to second order implies the usual Slavnov-Ward identities for the gravi- 
ton self-energy in the loop graph sector and generates the correct quartic 
graviton-matter interaction in the tree graph sector. The mass zero case is 
also discussed. 
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1 Introduction 



In this paper we follow the quantum field theoretical approach to gravitational 
interactions coupled to scalar matter fields (see the introduction to this subject 
in [Q] and reference therein). This approach allows a quantization of the in- 
volved fields, matter and graviton fields, and a Lorentz covariant perturbative 
expansion of the scattering matrix S. 

Calculations of matter loop diagrams in this conventional framework led to 
non-renormalizable ultraviolet (UV) divergences ||] . These were later confirmed 
by means of dimensional regularization and background field method, both in 
the massive || and in the massless ||], || case. 

The counterterms needed to cancel the divergences are not of the type 
present in the original Lagrangian density. According to these findings, quantum 
gravity (QG) coupled to matter fields does not fulfil the criterion of perturbative 
renormalizability ||. 

We show how it is possible to overcome these discouraging outcomes by 
applying an improved perturbation scheme which has as central objects the 
time-ordered products and as constructing principle causality. The S-matrix is 
constructed inductively as a sum of smeared operator-valued n-point distribu- 
tions avoiding UV divergences. 

This idea goes back to Stiickelberg, Bogoliubov and Shirkov and the pro- 
gram (causal perturbation theory) was carried out successfully by Epstein and 
Glaser 0] for scalar field theories and subsequently applied to QED by Scharf || , 
to non-Abelian gauge theories by Diitsch et al. ]|] and to quantum gravity [ 10 1 
in the last years. 

In addition, QG has considerable gauge properties ]ll] , which are formulated 
by means of a 'gauge charge' that generates infinitesimal gauge variations of the 
fundamental free quantum fields, Sec 2.3 . 

The present work focuses mainly on three aspects of QG coupled to massive 
matter fields. Brief remarks for the massless case are given for these aspects. The 
first aspect is the UV finiteness of loop graphs, which include the lowest order 
massive and massless scalar matter loop corrections to the graviton propagator, 
Sec. pi and the matter self-energy, Sec. ||. The result are UV finite and cutoff-free 
due to techniques of causal perturbation theory. 

The second aspect consists in the investigation of the gauge properties of the 
graviton self-energy, Sec. 2.4 . Gauge invariance of the S'-matrix implies some 
identities between the C-number parts of the n-point distributions which yield 
the gravitational Slavnov-Ward identities (SWI) fjffl , Sec. 3.3 . 

The third aspect of this work is also connected with gauge invariance: per- 
turbative gauge invariance to second order in the tree graph sector requires 
the introduction, at a purely quantum level, of a quartic matter-graviton in- 
teraction exactly as prescribed by the expansion of the classical matter-gravity 
Lagrangian, Sec. [B|. 

The quantization of the graviton field, the identification of the physical sub- 
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space and the proof of 5-matrix unitarity are investigated in |l3| , which provides 
also the conventions and the notations used here. Calculations involving gravi- 
ton self-couplings are not considered here, see [14|. 

The causal scheme applied to quantum gravity coupled to photon fields leads 
also to very satisfactory results with regard to the UV finiteness in loop calcu- 
lations and to the gauge invariance [15|. 

We use the unit convention: h = c = 1, Greek indices a,/3,... run from to 
3, whereas Latin indices i, j, . . . run from 1 to 3. 



2 Quantized Matter— Gravity System and Perturba- 
tive Gauge Invariance 

2.1 Inductive Construction of Two-Point Distributions in the 
S-Matrix Expansion 

In causal perturbation theory ]|], the ansatz for the 5-matrix as a power 
series in the coupling constant is central, namely S is considered as a sum of 
smeared operator-valued distributions 

°° 1 f 

S(g) = 1 + 5Z — i d 4x i---d A x n T n (xi,... ,x n )g(xi) ■ ... ■ g(x n ) , (2.1) 

n=l n ' J 

where the Schwartz test function g E 5(M 4 ) switches the interaction and provides 
a natural infrared cutoff. The S'-matrix maps the asymptotically incoming free 
fields on the outgoing ones and it is possible to express the T n 's by means of 
free fields without introducing interacting quantum fields. 

The n-point distribution T n is a well-defined 'renormalized' time-ordered 
product expressed in terms of Wick monomials of free fields : 0(x\, ... ,x n ): 

Tfi {x\ , . . . , X-n) — ^ '. 0{xi , . . . , Xn) . t n (yX\ X n , • • • , X n —\ Xn) . (2.2) 

o 

The tn are C-number distributions. T n is constructed inductively from the first 
order T\{x), which describes the interaction among the quantum fields, and 
from the lower orders Tj, j = 2, . . . , n — 1 by means of Poincare covariance and 
causality. The latter leads directly to a UV finite and cutoff-free distribution 
T n . 

For the purpose of this paper, we outline briefly the main steps in the con- 
struction of T^(x, y) from a given first order interaction. Following the inductive 
scheme, we first calculate the causal operator-valued distribution 

D 2 {x,y):=B! 2 {x,y)-A! 2 (x,y)= [T x {x), T^y)} . (2.3) 

In order to obtain D 2 (x, y), one has to carry out all possible contractions between 
the normally ordered T\ using Wick's lemma, so that D 2 (x, y) has the following 



4 



structure 



D 2 (x,y) = J2-0(x,y): d°(x-y). (2.4) 

o 

d^ix—y) is a numerical distribution that depends only on the relative coordinate 
x — y, because of translation invariance. 

D2(x,y) contains tree (one contraction), loop (two contractions) and vac- 
uum graph (three contractions) contributions. Due to the presence of normal 
ordering, tadpole diagrams do not show up. In this paper we do not consider 
vacuum graphs. D2(x,y) is causal, i.e. supper? (z)) C V + (z) U V~(z), with 
z := x — y. 

In order to obtain T2(x,y), we have to split D2(x,y) into a retarded part, 
R2(x,y), and an advanced part, A2(x,y), with respect to the coincident point 
z = 0, so that supp(i?2(-£)) ^ V + ( z ) an d supp(A2(z)) C V~{z). This splitting, 
or 'time-ordering', has to be carried out in the distributional sense so that the 
retarded and advanced part are well-defined and UV finite. 

The splitting affects only the numerical distribution d° (x — y) and must 
be accomplished according to the correct 'singular order' co^d^ ) which describes 
roughly speaking the behaviour of d^ix—y) near x—y = 0, or that of d% (p) in the 
limit p — > oo. If ijj < 0, then the splitting is trivial and agrees with the standard 
time-ordering. If u > 0, then the splitting is non-trivial and non-unique: 

d° 2 (x-y)—> rf ( x -y)+Y J C a , a D a «5< 4 > (x - y) . (2.5) 

|a|=0 

A retarded part r°{x — y) of d^ (x — y) is usually obtained in momentum space 
by means of a dispersion- like integral, see Eq. fl3.18|) . 

Eq. fl2.5|) contains a local ambiguity in the normalization: the C aj0 's are 
undetermined finite normalization constants, which multiply terms with point 
support D a 5 {i) {x — y) (D a is a partial differential operator). This freedom in 
the normalization has to be restricted by physical conditions. 

Finally, T2 is obtained by subtracting R' 2 (x,y) from R,2(x,y) and the whole 
local normalization coming from ( |2.5| ) is called N2(x,y). 

2.2 Quantized Matter-Gravity Interaction 

We consider the coupling between the quantized symmetric tensor field h^ v {x), 
the graviton, and the quantized scalar field (fi(x), the matter field, in the back- 
ground of a Minkowski space-time. 

The free scalar field of mass m satisfies the Klein-Gordon wave equation 

(□ + m 2 )(p(x) = 0, (2.6) 
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which follows from the free matter Lagrangian density Cfj = ^cf) }U (f>> 1 ' — ^-(fr 2 - 
The matter energy-momentum tensor reads 

tit : = w - tt c<$ = - -rr <m iP + < ^ > (2-7) 

where ij^ u = diag(+l, —1, —1, —1) and it fulfils Tj^ v = 0. Quantization of the 
scalar field is accomplished through 

[<t>{x),<j>{y)]=-iD m (x- y ), (2.8) 

where 

D m {x) = D£>(x) + £>£>(x) = Jd 4 p5(p 2 - m 2 ) sgn(p°) (2.9) 

is the causal Jordan-Pauli distribution of mass m. 
The free graviton field satisfies the wave equation 

nh^(x) = (2.10) 

and is quantized (see Q|) according to 

[h' w (x),h a ^y)]=-iV a ^D (x-y), (2.11) 

where the 6-tensor is constructed from the Minkowski metric 

and Dq(x) is the Jordan-Pauli distribution of Eq. ( ^S| ) with m = 0. 

The graviton field interacts with the conserved energy-momentum tensor of 
the matter fields. The first order matter coupling is chosen to be 

Tf(x)=i^:h^(x)b a p^T^(x):=i^{:h a ^ >a( t> > p: }, (2.13) 

where k is the coupling constant (see below for its relation to Newton's constant). 

To simplify the notation, the trace of the graviton field is written as h = /i 7 7 
and all Lorentz indices of the fields are written as superscripts whereas the 
derivatives acting on the fields are written as subscripts. All indices occurring 
twice are contracted by the Minkowski metric yf^ . We skip the space-time 
dependence if the meaning is clear. 

The presence of the 6-tensor in Eq. ( |2.13| ) is a consequence of the choice of 
the graviton variable (see Sec. 2.3). 
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2.3 Perturbative Gauge Invariance 

The classical gauge properties of h^ u {x) (which are related to the general co- 
variance of the metric g^ u {x) under coordinate transformations [17], [|l^]) are 
formulated at the quantum level by the gauge charge [10], [11| 



Q:= J &xhr{x)p&Z u^x). (2.14) 

x°=const 

For the construction of the physical subspace and in order to prove unitarity of 
the ^-matrix on the physical subspace [l3|| , the ghost field (x) , together with 
the anti-ghost field u v {x), have to be quantized as free fermionic vector fields: 

□«'*(z) = 0, Du v (x)=0, {u»{x),u v (y)} = irf v D (x - y) , (2.15) 

whereas all other anti-commutators vanish. The gauge charge generates the 
following infinitesimal operator gauge variations 

dqhTix) := [Q,h^(x)] = -ib^uP{x). a , 
d Q u a (x) := {Q,u a (x)} =0, d Q u a (x):={Q,u a (x)}=ih af} (x) t p. (2.16) 

Gauge invariance of the 5-matrix means formally 

lim d Q S(g)=0. (2.17) 

This condition can be reformulated in terms on the ?vpoint distributions T n \ 
using Eq. ( j2.1|) , we see that the condition of perturbative gauge invariance to 
n-th order 

dgT n (xi, . . . , x n ) = sum of divergences , (2-18) 

implies Eq. ( |2.17| ), because divergences do not contribute in the adiabatic limit 
g — > 1 due to partial integration and Gauss' theorem. 

Using a simplified notation which keeps track of the field type only, per- 
turbative invariance to first order in pure QG |uj, namely for a coupling of 
the form T± ~: hhh : without matter fields, requires the introduction of the 
ghost coupling T" ~: uhu : , so that d Q {T[ l + T?)(x) = d^T^ +u (x). Here, 

Tyi +U ^:uhh: + :uuu: is the so-called Q-vertex. 

Using Eq. (2.16), the first order matter coupling (2.13) is gauge invariant: 

dQT?(x) = \b^° :uP{x), tJ b a ^ v T^{x):= \ :u"{x), a T^(x): 
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{ :u<>^y. -~ :vTMy. :u^<j>: }) ( 2 - 19 ) 

r I 
A 



-■d§T^(x), 
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because of T^ a = 0. Ty^(x) is the matter Q- vertex. The concept of Q- vertex 
allows us to formulate in a precise way the condition fl2.1S|) of perturbative gauge 
invar iance to the n-th order: 

n d 

d Q T n (xi,... ,x n ) =^2 — T^ /l (x 1 ,... ,xi,... ,x n ), (2.20) 
1=1 %l 

where T^n is the 'renormalized' time-ordered product, obtained according to the 
inductive causal scheme, with one Q-vertex at X[, while all other n — 1 vertices 
are ordinary T\ -vertices. 

The procedure outlined here corresponds to the expansion of the Hilbert- 
Einstein and matter Lagrangian density 0, ||, || 

Ceh + £ M = ^ V=g 9^ R»v + Iv^g {aT^v - ™ 2 4?) , (2.21) 

(k 2 = 32-7T G) written in terms of the Goldberg variable g^ u , in powers of the 
coupling constant k according to the metric decomposition 

$>» : = y/^jjg"" = rf v + k h? v , (2.22) 

which defines the graviton field W in the Minkowski background. Then one 
obtains 

00 

C-EH + C'M = ^ {/^E!H ~l~ ) • (2.23) 

3=0 

From choosing the Hilbert gauge = 0, one obtains Eq. ( |2.10| ) and the 
presence of the 6-tensor is made clear fLsfl . 

The first order graviton coupling T^(x) ~: hhh :, corresponds then to the 
normally ordered product of ikjCbh (see for a derivation based merely on 
the principle of perturbative operator gauge invariance) and C^h ~ hhhh rep- 
resents the quartic graviton coupling required by perturbative gauge invariance 



to second order in the tree graph sector [10| 



The expansion of the matter Lagrangian density reads 
Cm = \ (tTMv ~ m 2 cf> 2 ) + ~ - ^ %i/ </> 2 ) 



-(0) - K r { - 1 ^ 

'M M 



+ V^L (h^h^M - I hh^cp) +0(k 3 ) . (2.24) 
8 2 



From the first term one obtains (2J3) and the matter coupling of Eq. ( 2.131 ) 
corresponds to ik : :. A quantized quartic interaction ~: hh(fxf) : which 
agrees with C!\} will be necessary for reason of gauge invariance, see Sec. [|. 
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2.4 Identities for the Two-Point Functions from Perturbative 
Gauge Invariance to Second Order 

From the structure of it is evident that the two-point distribution describing 
loop graphs has the form (up to non-contributing divergences for the matter self- 
energy, see Sec. 4.1): 

T 2 (x,y) l °°v s =:h a P(x)h^(y): iU(x - y) a ^ v + :<t>(x)ct>(y) : iE(x-y). (2.25) 

Here, the first term represents the matter loop graviton self-energy and the 
second term the scalar matter self-energy. The C-number distributions U(x — 



y)aj3jj,u and E(x — y) will be explicitly calculated in Sec. ^2 and in Sec. [4.2| , 
respectively. 

Perturbative gauge invariance to second order, namely Eq. (|2.20| ) with n = 2, 
allows us to derive a set of identities for these numerical distributions by compar- 
ing distributions attached to the same operators on both sides of Eq. ( 2.20D , [19|. 

We compute d Q T 2 (x, y) loo P s by means of ( |2,16D and isolate the contributions 
with external operator of the type :u(x)h(y) :. We obtain 

dQT 2 {x,y)^% u{x)h{y) . =:uO{x), a h^{y): (b a ^ U(x - y) a ^ u ) . (2.26) 

On the other side, T^ ± (x,y) has to be constructed with one Q-vertex at x and 
one 'normal' vertex at y. From the structure of both interaction terms, it follows 
that the loop contributions coming from T^nix^y) can only be of the form 



L 2 



^ /1 (x,y)=:uP(x)h^(y): ^(x - y) p ^+ :u°(x)h^(y) : t uh (x-y)^, (2.27) 



by performing two matter field contractions. The second term T^ 2 (x,y) does 
not contain terms with Wick monomials of the type : u(x)h(y) :. Applying d% 
to the expression above we find 

d x a T° /1 (x,y)=+:u»(x), (7 h^(y): {t^(x - y) mu + r?/ t uh (x - y )^} + 
+ ; U P{x)h^{y): dl{e uh {x-y) p ^ + r ] °t uh {x-y), u }. 

We compare the C-number distributions in ( |2.26| ) and in ( |2.2§| ) attached to 
the external operators 

:u p {x)^ u {y): and :u p (x)h^(y): . (2.29) 

Therefore, we obtain two identities 

b?™? n(x - y) a ^ = {t° uh (x - y)% + rT t uh (x - y), u ) , 

= d*{t° h (x-yy^ + r l <>°t uh (x-y)^}. 

By applying to the first identity and inserting the second one, we obtain 

b a ^° d*Tl(x - y) a ^ u = . (2.31) 



9 



This identity for the matter loop graviton self-energy tensor has been explicitly 
checked and implies the gravitational Slavnov-Ward identities for the two-point 
connected Green function (see Sec. 3,3). 

Gauge invariance to second order in the tree graph sector is much more 
involved and requires also the full analysis of the matter-graviton interaction, 
Sec. |. 
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3 Matter Loop Graviton Self-Energy 
3.1 Causal D 2 (x, ^-Distribution 

In order to construct D2(x,y), according to Sec. we first need the contrac- 
tions between field operators. From fl2,§| ) and ( |2.11| ), we derive them: 

C{</>(x)</>(y)} := [<l>(xY-\<l,(y)^]=-iD^\x-y), 
C{h a/3 (x) h^(y)} := [h ap {x) ( -\ h^(y) (+) ] = -i b a ^ u D ( +) (x - y) , (3.1) 

where (±) refers to the positive/negative frequency part of the corresponding 
quantity. 

The A' 2 (x,y) 9SE distribution for the graviton self-energy by a matter loop 
is obtained by performing two matter field contractions in — T 1 M (x)T 1 M (y). Us- 



ing (3.1) and with = —da-, we find 

.2 



A f 2 (x,yr E = :h^(x)h^(y): ~ A 



FP r) x n<+> . rFrP /")(+) 4- 



4 L 

- &*J? V D<$ ■ SFffiD™ + m 2 V, d-DM . 8* D W + (3.2) 

4 

+ m 2 VaP d?D%> ■ 8%D% - ^- riopn^ • DW] (x-y). 
We introduce the functions 

D^ m (x) := D%\x) • DW(x) , D^f(x) := 8%D^(x) • ^(j) , 

■= KtfDgix) ■ &Z82Dg{x) , (3.3) 

so that we have 



A' 2 (x,yy SE =:h^(x)h^(y): a' 2 (x - y) 9 a s E 



L> (±) ,™ - D 



(±).m | 2 T-\fi),m i 

Jw + m V^ D a\p + 



I 2 T^(±),m m T->(±),m 



(3.4) 



(x-y). 



Products of Jordan-Pauli distributions are evaluated in momentum space with 
Eq. ([2.9|), because there products go over into convolutions between the Fourier 
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transforms. For example, D\'' m (x) becomes 
-1 



£>y ),m (p) = —^4 jd 4 k5((p - kf - to 2 ) 6( ± (p° - A; )) ,5(£; 2 - m 2 ) G(±£; ) . 

(3.5) 

Therefore, the basic integrals that remain to be computed are of the form 
I^(p)~/a/afS/a^/a^ := fd 4 k 5((p - k) 2 - to 2 ) Q( ±(p° - k )) 5(k 2 - to 2 ) 



x 0(i/c ) [l, k a , k a kp, k a kpk^, kakpk^k^ , 



(3.6) 



which are calculated in App. 1. By means of the Im (p)... -integrals, the D , - 
functions in momentum space are 



DT hm (p) 
+1 



1 



(2tt)< 



D ( ^f (P) = 7^4 [ + Pa I™ ( P )f3 ~ 4 ±5 (P W 



D 



(±) - m ( v ) = L r 



+ PaPf) I { m {p)nv ~ Pa Iffi (p)/3^ 
PP J m (P)a^ + (p)af3»v ■ 



(3.7) 



Inserting Q with Eqs. flO|) , flAj), flOU) , ( [Alol) and (gig) into Q, then 
the a^-distribution in momentum space reads 



-K 2 TT 



" 960(2vr) 4 



aPpPfiPv 

+ B p 2 {paPpripv + p^Pullap) + 
+ C p 2 (paP^VPv + PaPuVfSfi + PpP^Vav + PpPuVafi) + 

+ £ p A {q afl -r]i3u + rjca/V/Sit) 

+ Fp 4 Vaf3V f iv dip)™ , 



(3.8) 



with the coefficients 



2 4 
TO , TO, 

A := -8 - 16 - 48 — r 



B :=-4-8 



TO 



24 ■ 



P 



TO 
4 

r 



p 

P 



2 4 2 

„ TO TO „ TO TO 

C:= +1-8 -=- + 16-3-, E := -1 + 8-^- - 16 — 
F:=-l-12^+4^, 



(3.9) 
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and the distribution d(p)^ := w 1 — ^p- 0(p 2 — 4m 2 ) 0(±p°). Performing the 
same calculations for R' 2 (x,y) = — T 1 M (y)T 1 M (x), we obtain 



-K 2 7T 



^(p)Si/ = 96o( 27r )4 t the same as in Eq - ©] d (P)™ ' 
Therefore, with (|2.3|) the causal D 2 (x, y)-distribution reads 



D 2 (x,yr E =:h^(x)hi"(y): d 2 (x - yffi 



K 2 7T 



960(2vr) z 



[the same as in Eq. Q3.8] )] d(p) r 



Here, d(p) m = d(p)« - d(p)^ = ^1 
distribution can be recast into the form 



4m 2 



(3.10) 



(3.11) 



6(p 2 - 4m 2 )sgn(p°). The d 2 - 



i=l 




2 4 
777- * ?7t ~ 1 

P(p)a/3tiu + Q(p)a/3fiu + ~^ R(p)a/3fiv d(p) r 



(3.12) 



where the polynomials of degree four are given by their coefficients 

P(p)a^u = [-8, -4, +1,-1,-1] , 

Q(p)a/3^ = [ ~ 16, -8, -8, +8, -12] , 
R(p)a^u = [ ~ 48, -24, +16, -16, +4] , 



(3.13) 



according to the structure given in Eq. 

In the case of massless (m = 0) matter coupling, that is the first order 

matter interaction is chosen to be T™(x) = : h a/3 (x)(f)(x) >a (j)(x) ^ :, then the 
^-distribution reads 

4(p)^=° = T P(p) a ^ &(p 2 ) sgn(po) • (3.14) 

Hence, the limit m — > of ( |3.12 ) is feasible without problems, see Eq. ( |3.48| ) for 
the splitting in the m = case. 

The extension to non-minimally coupled massless matter is also considered. 
From 



(3.15) 
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K ( 2 

2l 3 



we derive the first order matter coupling 

which yields (see [|14|], |15|] for the calculations in the m = case) 
rf 2 (p)^™' =X[-H,_| j i, -l, «] (p 2 ) sgn(po) • 



The corresponding i2-distribution will be given in Eq. (|3.49|). 



(3.16) 



(3.17) 



3.2 Distribution Splitting and UV Finiteness for the Matter 
Loop 

We turn now to the splitting of the /^-distribution of Eq. ( |3.12| ). The leading 
singular order is four because the polynomials are of degree four in p. But, 
since these polynomials act in configuration space as derivatives, the essential 
structure of the distributions is given by the scalar part. Therefore, neglecting 
the polynomials, the first, the second and the third term in ( 3.12| ) has singular 
order 0, —2 and —4, respectively, due to the inverse powers of p. When discussing 
the normalization N2(x,y), we will realize that this was the correct choice. 

As anticipated in Sec. ^]l], a retarded part of c?2 is obtained in momentum 
space by the integral M: 



dkn 



d(k ) 



(ko-iO^ipo-ko + iO) 



(3.18) 



for p v = (po,0), po > 0. This retarded part is the so-called 'central splitting 
solution', because the subtraction point is the origin. 

The behaviour of the first term in ( |3.12|) is dictated by a scalar distribution 
of the form 



d(k) 



(i) 



1 



4m 2 



0(k 2 - Am 2 ) sgn(fco) =: f(k 2 ) sgn(A; ) . 



(3.19) 



Therefore, we have to split d(k) with uj(d) = 0, because d(ko) ~ 
\ko\ -» oo and k v = (k ,0) £ V+. From ( |3l8| ) we obtain 



constant for 



i 




r+oo 




PoJ 


1 dk 


~ 2^ 


— oo 


i 




POO 


~ 2^ 


Po 


/ dk 

Jo 


i 




POO 


~ 2^ 


POJ 


1 dk 




i 




POO 


~ 2^ 


POJ 


1 dk 





/(fco)sgn(fc 



(fc - «0)(po - k + tO) 
f(k 2 o) 



ko(Po ~ k + i0) 
f(k 2 ) / 1 



f(k 2 ) 



dko 7 

oo k {po - k + t0) 



(3.20) 



+ 



ko ^Pa - ko + i0 po + ko + iO 
fig) 2pok 
k o Po~ k o + *Po° ' 
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With the new variable s := fcg we find 



r c (po) = f pg f& , s /( *? . m • (3.21) 



Inserting the explicit form of f(s) we have 



r c (Po) = ^-Po I ds\ll — 2 — — -. (3.22) 

2vr J 4m 2 V s s(Pq - s + zpoO) 

We decompose the integral into real and imaginary part according to (x+iO)^ 1 = 
P(x~ l ) -inS(x): 



r c (po) = ^PoPl ds \l l -———nr — \ + 



4m 2 

2tt ^ u ~ J Ami" V ~ s s(pl - s) 



1 / /J. ?77 2 

+ 2 A/ 1 --^-©(Po-4m 2 )sgn(p ). (3.23) 

The T2(x, y)-distribution is obtained from the retarded part B,2(x,y) by sub- 
tracting R' 2 (x,y). This subtraction affects only the numerical distributions. 
Therefore, subtracting 



f'(po) = -Jl ~ ^ e(pg _ 4m 2 ) e(-po) , (3.24) 
V Po 

from ( |3.23| ), we obtain the numerical f-distribution belonging to T2(x,y): 
i{Po) = r c (po) ~ r{po) 



oc 



4m. 2 1 1 / 4m 2 

(3.25) 



-piP I ds\ll-^——^ - + - x ll- ^-Qipl- Am 2 

2tt P ° 7 4m2 V a s(p 2 -s) 2 V p 2 ^° 



which can be written in the form 

i 



f°° / 4m 2 1 

J J-r—r , w -.+«) • <3 ' 26) 

This result can be generalized for p 6 y + by introducing the 'inverse momentum' 
q := 4m 2 /p 2 so that 



Note that we write for simplicity the p-dependence instead of the g-dependence 
of the basic integral II that remains to be calculated. Therefore, the splitting of 
the first term in ( 3.12j ) and the subtraction of yields 

Kp) { X„ = * 2 HpU^u fl(p 2 ) ■ (3.28) 
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Here, E := T/(2vr) = K 2 vr/(960(27r) 5 ). 

Following the same steps as from Eq. ( 3.2C| ) to Eq. ( |3,28| ), we find for the 
second term in ( |3.12j ) 

2 

= * S Q(pUv» ft (p 2 ) , (3.29) 
because in d(p)^L = m? T Q{p) a pia/ d{pY 2 \ the scalar part reads 



Analogously, for the third term of ( 3.12j ) we obtain 



777 

Kp)%, v = ^~^ R(pU,u J(P 2 ) , (3.31) 



because in d{p)^p = m 4 T R(jp) a piJ,v d(pY 3 \ the scalar part reads 



d(p) {3) = 4 \h - ^ ®(P 2 - 4 ™ 2 ) sgn(po) , (3.32) 



pA W p2 



with the definition 



J(p 2 ) := f°cto s ^ (S ~ g) n T • (3.33) 
V 7 7 g s 3 (l-s + z0) v ; 

The difference between J(p 2 ) and fl(p 2 ) lies in the powers of s in the denomina- 
tors in Eqs. (^27|) and ( fOp . This is a consequence of splitting these distribu- 
tions according to the singular orders of the corresponding scalar distributions. 

The two integrals tl(p 2 ) and J(p 2 ) have the same structure and the last can 
be expressed by means of the first. We decompose Tl(p 2 ) into real and imaginary 
part 

tl(p 2 ) = P f°°ds V 2 S [ S ~ g) - i vr 71^6(1 - q) , (3.34) 

and concentrate our attention to the principal value part. With the substitution^] 
s(s — q) = (s — x) 2 the real part of II(p 2 ) reads 



/•oo 

f[ r (p 2 ) = -2P dx 

Jo 



(x - qf 



x 2 (x 2 -2x + q) 

op °° d J± 1-g ^ ( °- >0) 

x 2 x 2 — 2x + 



1 we choose a;(s) = s + sj s(s — q), so that x(s) goes from q to oo for s going also from q to 

oo. 
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having factorized the integrand. The j(p 2 )-integral yields a real part 

(x - q) 2 (2x - q) 



Jr(p 2 



f-OO 

-2P I dx 



2P dx 



x 4 (x 2 - 2x + , 



1 2g_ f_ 



1 



x s x 4 x 2 — 2x + q 



(3.36) 



A look at Eqs. (3.35) and ( |3.36|) en ables us to isolate in the expression for j r (p 2 ) 
the terms appearing also in ( |3.35| ). The others can be easily integrated and we 
obtain 



Mp 2 



3q 



+ n r ( P 2 



6m 2 



+ n r ( P 2 ), 



(3.37) 



The imaginary parts have always the same form as in Eq. ( 3.34 ). Gathering 
the results in ( 3.28| ), ( |3.29| ) and ( 3.31| ) with ( |3.37| ) we can write the distribution 
describing the matter loop graviton self-energy: 



cefifiv 



i=l 



{Pip) 



V 



P 



• <»" + 5- Q(p)<*piw + 3r RiPaPi**)} n(p 2 ) + ( 3 - 38 ) 



+ 



m 
Qp' 



R{p)af3fiv 



Therefore, the 2-point operator- valued distribution T%(x, y) for the graviton self- 
energy reads 



gSE r a/3 



h^{x)h^{y): t 2 {x-y)%l u =:h^{x)h^{y): iU(x - y) a ^ , 

(3.39) 



where Il(x — y) a pu,v is the graviton self-energy tensor. Its Fourier representation 
is given by — zx ( |3.38| ). 

We still must calculate explicitly the integral representation for II (p 2 ), (3.27). 
There are three different regimes, depending on the value of q. For q = 1, namely 
p 2 = 4m 2 , we obtain by means of the partial decomposition (|3.35|): 



fl a (p 2 = 4m 2 



OO 1 

dx — =■ 



-2. 



(3.40) 



For q < 1, namely p > 4m , the integration of the partial decomposition (|3.35| ), 
taking into account also the imaginary part from ( 3.34 ), yields 

fl fe (p 2 ) = -2 + ^T~q lo. 



q-l-VT 



q-l + ^T^q 



ny/T=qe(l-q). (3.41) 
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For q > 1, namely < p 2 < Am 2 , the integration of (3.34) gives 
fl c (p 2 ) = -2 + 2 sfq~=l 



7T q- 1 

arctan , 

2 



= —2 + 2 yjq—l arctan 
Note that these three results are connected by 



1 



(3.42) 



lim U b (p 2 ) = U a (p 2 = Am 2 ) , lim Ii c (p 2 ) = Ii a {p 2 = Am 2 ) . (3.43) 

p 2 \4m 2 p 2 /4m 2 

Writing the p-dependence explicitly, the final form for fl(p 2 ) is 

1 - 



+ W1 



Am s 
p 2 



log 



1 



4m 2 



1 



1 



4m 2 



+ i 7T 



-2 



/ 4 m 2 



1 arctan 



4rrt 2 



Q(p 2 -Am 2 ) + 
@(Am 2 -p 2 ) I . 



1 



(3.44) 



Two limits of this result will be used in the discussion of the normalization 
N2 of the T2-distribution, Sec. |3.4| : the limit of tl(p 2 ) for p 2 \ and the limit 
of tl(p 2 )/p 2 for p 2 \ 0, too. In the first case we have 

lim fl(p 2 ) = lim Tl c (p 2 ) 

p 2 \0 q-^oc 



lim 

q— too 

lim - 

q— >oo 



2 + 2Jq~^l (--(-- + + 0(- 



,5/2- 



2 + 2 



3(9 - I! 



+ 0( 



For the second limit, we obtain 



p 2 \0 



lim 



(3.45) 



(3.46) 



lim - — ? 

g— >oc 4m z 



{ _ 2 + 2v ^T (_J_ - 4- Ot-lj)) } = ^ . 

The last consideration concerns the retarded part in Eq. ( p. 23 ), given also 
by ( fT44| ) up to the signum-function in j»o : this retarded part is the boundary 
value of the analytic function 



r(p) c 



Am 2 
p 2 



log 



^ 4m 2 ^ 



^ 4m 2 1 ^ 



(3.47) 
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Summing up the whole calculation, we have found the 2-point distribu- 
tion ( |3.39D for the graviton self-energy contribution. The corresponding ten- 
sor, the structure given by (3.38) and the integral in (3.44), is UV finite and 



cutoff-free. 

During the calculation, we never resorted to an ad-hoc regularization of the 
expressions (for example dimensional regularization as in 0). This was made 
possible by using the correct starting point, namely Eq. ( |3.18| ), which is, so to 
say, a careful multiplication by a step- function in the time argument. If this 
had been done naively, then it would have corresponded to the choice uj = — 1 
in ( 3.18j ), when splitting the first term of Eq. ( |3.12 ), a choice which is manifestly 



wrong, being uj = 0. 

Choosing uj = —1 in Eq. (3.20), one obtains a UV logarithmic divergence. 

For the sake of completeness, we briefly report also the results in the case of 
massless matter coupling, Eq. ( p. 14 ), and in the case of non-minimally coupled 



massless matter, Eq. ( 3.171) . 

The splitting of the scalar distribution 0(p 2 ) sgn(po) requires some modifi- 
cations if one tries to use the splitting formula ( |3.18[ ), see Q, [14]. The retarded 
part is given by (i/27r)log ((— p 2 — ipoO)/M 2 ^), so that the m = matter self- 
energy tensor reads 

U(p)^=l = E P(pU, v log (^j^) , (3-48) 

where M > is a scale invariance breaking normalization constant and not a 
cutoff. 

For the non-minimally coupled case we find analogously 

n(.)~ = S [- H, X) -1, «] log (=^2) . (3.49) 

This graviton self-energy tensor is traceless: rf 1 ^ ^(p)^ 1 ^™ 111 ' = 0, because in this 
case the graviton is coupled to a traceless matter energy- momentum tensor. The 
latter corresponds to the so-called 'improved' energy-momentum tensor |2(J. 
In addition, it is transversal: p a n(p)„^™ in ' = 0. This property follows from 

the gauge identity (|2.31|) , namely b 01 ^^ p a tl{p) a p llv = 0, see Sec. 3J3, and from 
its vanishing trace. 

Also in these two cases, we have found UV finite graviton self-energy contri- 
butions without introducing counterterms or UV cutoff. This is in contrast to 
the calculations carried out for massless scalar matter fields coupled to QG in 
the background field method with dimensional regularization 0], ||. 
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3.3 Graviton Self-Energy Tensor and Perturbative Gauge In- 
variance 

The gauge properties of T2(x,y) 9SE are contained in the identity b a ^ pa d§H(x — 
y)ai3^v = 0, see Eq. ( 2.31 ). This identity implies the conditions 

71-25 = 0, C + E = 0, B -2E -2F = 0, (3.50) 

for the coefficients of the self-energy tensor. These conditions are satisfied by 
our result of Eq. (|3.9[) and therefore II(x — y) a /3fj,u is gauge invariant. This is 
certainly the case at the level of the D2(x, y) 5SB -distribution, before distribu- 
tion splitting. In the causal scheme, UV finiteness and gauge invariance are 
established separately. The latter is not used to reach the former. The iden- 
tity ( 2.31 ) implies the Slavnov-Ward identities (SWI) for the 2-point connected 
Green function. The latter is defined as 

G{p) [ Xu ■= W AT (p) fl( P y 5pCT W D F Q (p) , (3.51) 

where D F (p) = (27r) _2 (— p 2 — iO)^ 1 is the scalar Feynman propagator. The 
two attached lines represent free graviton Feynman propagators. The SWI 
reads §, @: 

P a G(p) [ Xu = > (3-52) 

namely that the 2-point connected Green function is transversal. In term of the 
coefficients A, . . . , F as in Eq. ( |3.§| ) we have 

A-2B = 0, C + E = 0, - + C-F = 0. (3.53) 

These are equivalent to the conditions ( |3,50| ). This conclusion is valid also in 
the massless case, Eqs. (|3l8|) and (gl|). 

In [^TJ , the gauge invariance of the massless matter loop graviton self-energy 
tensor is also investigated, but there it is not realized that the correct matter 
coupling is the one in Eq. ( |2.13 ), namely with the 6-tensor, when one uses the 
Goldberg variable expansion. This deficiency does not have consequences, if 
the graviton is coupled to a traceless matter energy-momentum tensor as in the 
non- minimal coupling case, Eq. ( 3. 16| ) . 

The condition of perturbative gauge invariance to second order in the loop 
graph sector 

d Q T 2 (x, y y SE = S£( :uP{x)h^(y) : [b a ^U(x - y) a ^ v ]) + (x <-> y) 
= d x a ( K :uP{x)h^{y): [^(x - y) pia/ + ^ t uh (x - y)^] ) + (x <- y) , (3.54) 

has been explicitly checked by calculating also the distributions t° h (x — y) Pilv 
and t u h(x — y)^ with one Q- vertex from Eq. ( |2.27 ). 
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3.4 Reduction of the Freedom in the Normalization 



We turn now to the normalization of the T2-distribution of Eq. ( 3.3S| ). The total 



singular order is four because the polynomials are of degree four in p. Therefore, 
we have to add normalization terms up to the singular order four. 

The freedom in the normalization due to the splitting procedure is contained 
in the local term N 2 (x, y) 9SE : 

N 2 (x, y y SE =:h^(x)h^(y): i N(d x , d y ) a p^ 5™{x - y) . (3.55) 

From (|2.5|), we can write in momentum space this normalization as a sum of 
polynomials of degree 2i < 4 = to 

2 

N(p)a^=Y,^(P)al,- (3-56) 
i=0 

Normalization terms with odd uo are absent due to parity and Lorentz covariance. 
Gauge invariance b a ^ pa p lT N{p)^ = 0, (i = 0, 1, 2) and symmetries reduce the 
freedom in the normalization in such a way that the polynomials have to be of 
the form 

N(P) ( X„ = , N(p)%^ = E [0, 0, -a, a, -a] 1 , 

NipfXv = 3 [4(6 + c), 2(6 + c), -b, b, c] , (3.57) 

in the usual representation given by Eq. ( |3.8[) . The constants a,b,c € M 
should be fixed by requiring the appropriate mass- and coupling constant- 
normalizations for the corrections of order k 2 to the graviton propagator. Let- 
ting formally g — » 1 in Eq. (p.l|), we write the order k 2 corrected propagator 
as 



-iD{x - yYap^ = -ib al 3 f w Dq(x -y) + Jd i xid 4 x 2 ( -ib a p pa Dq(x - xij) 
x i (n(xi - x 2 y^ & + N( Xl - x 2 r^ 5 ) ( - i b jSfiV D F (x 2 - y)) . (3.58) 
In momentum space, this becomes 



= : n(p) a /3 M „ 

(3.59) 



After a little work, we find in the form of Eq. (| 

n(p)^ = h [A(p 2 ), / S ( P 2 ), / c (p 2 ), f E ( P 2 ), Mp 2 )} , (3.60) 
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with 

f A (p 2 ) 
fsip 2 ) 

Up 2 ) 



ar . , in ' \ 2 48m 



,.„ 2 m 4 \ 

16— -48— U(p 
p* p 4 / 

+6 + 32 — + 16 -J- n(p 2 ) + 



6p2 +^ + c) 



P 



P 



6p 2 



-ZAP 2 ) , 



+ i-8^. + i6^) n(p 2 ) + 



773 



16m 2 
6p 2 



With the formula (if g(p 
1 1 



^ - K 2 1 



1 



1 



2a 

26 - 4c + ^ 



6-^ = -/,b 2 ). 



1 



(3.61) 



> ' n ' > - n 9i.P 2 ) 2""^ = 2 ^ T^T +0 ^ 4 ) 

— iU —p A — i(j — |r — iU — — i\j —p A — i(j — gyp ) 



(3.62) 



we obtain the order k 2 corrected graviton propagator 



D{p) 



[2] 

aftfiv 



1 1 



+ ... 



(3.63) 



where non-contributing terms between conserved matter energy-momentum ten- 
sors have been neglected. The corrected propagator above has the correct limit 
lim^o-DO)^ = bap^Dfiip). 

Mass normalization (the graviton mass remains zero under quantum correc- 
tions) yields 



p A f E {p 2 



p 2 =0 



o = p 4 Up 2 ) 



p 2 =0 



(3.64) 



Since fl(p 2 = 0) = from Eq. fl3.45|) , these conditions always hold. 

Coupling constant normalization (k is not shifted by the quantum correc- 
tions) implies 



p 2 Up 2 ) 

Analysis of the first condition 



p 2 f E ( P 2 ) 



p 2 =0 



p 2 =0 



-16m 



o = p 2 Up 2 



p 2 =0 



U(p 2 



P 



16m 



+ a = 0. 



(3.65) 



(3.66) 



p 2 =0 



-l/6m 2 
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yields a = 0. Analysis of the second condition 



p 2 f F (p 2 



-16m 




16m 



2a = 0, 



(3.67) 



yields also a = 0. Decisive is the compensation between the first two terms 
in ( |3.66D and ( 3.67 ), This is due to the presence of the term ^iR{p)af3fjLv in 
Eq. (p|). 

A remark about the splitting: if we had split the distributions d^(jp)^Q , 
i = 2, 3 in Eq. ( |3.12| ) according to their 'true' singular orders, namely 2 and 
(because of the presence of the polynomials), respectively, then the term 
^R(p)a/3ixu would have been missing from ( |3.38 ). Working out the conse- 



quences for what concerns the normalization question, Eq. ( 3.66 ) would have 
required the choice a = — 8m 2 /3, whereas Eq. ([3.67 ) the choice a = 4m 2 /3. This 
would have meant the impossibility of a consistent normalization. Therefore, 



the splitting, as carried out in Sec. |3.2| , is justified. 

The origin of the above mentioned problem lies in the fact that the central 
splitting solution ( 3.18] ) is not applicable in that case and one has to choose a 
subtraction point different from the origin. 

The remaining constants b and c are not fixed by these requirements. The 
total graviton self-energy tensor including its normalization has then the form 



2 A 
TYl * TTt ~ 1 ^ 

P{p)a^v + Q(p)af3ixv + ~ R(Pa(3^)\ U (p 2 



P 



P 



+ 



m 



6p 2 



+ — R(p) a ^v + Yl ZiZ^afov 
i=l 



(3.68) 



where ZjGl,t=l,2 are still undetermined constants. The ^{pfap^v & are basis 
elements in the two-dimensional space of gauge invariant polynomials of degree 

\(2) 



a,uge invariant polyn< 
four. They can be chosen to be: Z{p) 1 ^ = [4,2,-1,1,0] and Z{p)f 



[4,2,0,0,1]. 

Analysis of the issue of normalization with the method used in [14|, [15| 
leads to the same conclusions. 



4 Matter Self-Energy 

The 2-point distribution describing the matter self-energy graph is not interest- 
ing from the point of view of its gauge properties. However, the calculation of 
the corresponding distribution is carried out to show its UV finiteness. 
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4.1 Causal D 2 (x, ^-Distribution and Distribution Splitting 

The ^-distribution is here obtained by performing one matter field contraction 
and one graviton contraction (|3.1|). The result reads: 



D 2 (x, y y 



--+ 

+ :<p{x 

+ :4>{x 
+ 



:<f>(x), a 4>(y), a : { (C<f 5 - C^) (x - y)] =: d a (x - y)} 



2 

-2^,2 



_2^2 



My),*'- {[A^{C ( .^-C^)(x-y)\ =:d b (x-yr} 

:d c {x-y) a ] 



urn-- {[^{&*>-c$){z- y j 

■4{x)ct>(y) :{[- K 2 ro 4 (C<+> - C^)(x - y) 

where the C (±) -functions are defined by 

C.f(x) :=D^\x)-D%\x), 
C${z) :=DM{x). 3K»(x). 



(4.1) 



(4.2) 



These products are calculated in momentum space, see App. 2, so that the 
distributions di, i = a,b,c,d read 



da(p) 



1 2 

-K m ir 



4(2vr) /i 



m 



d b (p) a 



• 2 2 

-in m w 



8(2vr) 4 
ddip) = 



m 



2 4 
K m IT 



2(2vr) 4 



1 



1 j- ) p a 0(p - m ) sgn(po) 



m 
p 2 



0(p 2 — m 2 ) sgn(po) • 



(4.3) 



From power-counting arguments, one could expect that d a behaves as p 2 for large 
momenta. This is not the case, because the wave equation (□ + m 2 )Dm\x) = 
lowers the power of p coming from the product of contractions. In order to 
shorten the calculation, we bring D2(x,y) mSE into the form 



D2(x,y) 



lllHE 



-: 4>(x)4>(y) : d 2 (x — y) m E + divergences , 



(4.4) 



with 



d2( P y 



+ P 2 d a {p) - ip a d h (p) a + ip a d c {p) a + d d {p) 




3m 2 m 4 



2 + 2p 2 



Q(p 2 — m 2 ) sgn(po) 



(4.5) 



Truly, this simplification can only be made for the corresponding T2-distribution, 
because divergences do not contribute in the adiabatic limit g — > 1 of Eq. (^l|). 
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Therefore, one should split the four distributions di, i = a,b,c,d separately and 
then recast the tj's in a form similar to Eq. (4^) for the T2-distribution. The 
final result would be the same. Note that in the m = case, D2(x,y) mSE = 0. 

The splitting of ( |4.5[) is accomplished by means of the splitting formula ( 3.18| ) 
with uj(d% SE ) = 2: 

i f' + °° Hn(kn^ mSB 



r c ( P o) mSE = ^ V / dko jr ," Q: u; - , (4.6) 

for = (po, 0), po > 0. The retarded part then reads 



»T , f + °° 17 9(Jfeg -m 2 )sgn(fc n ) L 3m 2 m 4 
2vr y.oo fc (p 

With s = fcg, ds = 2kodko, we obtain 



f r _ imS E_ » a r 00 ,. 0(^0 -m 2 )sgn(fco) f ovn „. , 



„ (>1 if 4 f°° 6(s-m 2 ) f 1 3m 2 ///' 



W" - ^ I * (p g! s + iw0) t; ""ST + 2?> • (4 ' 8> 

which can be decomposed into real and imaginary part: 



+ 2 e (Po - ™ 2 ) sgn(Po) |po - — + ^2 j • ( 4 - 9 ) 
Subtracting the distribution 



r' 2 (Po) 



o 3m 2 m 4 1 



e(p§-m 2 )e(- P0 ), (4.io) 



(coming from R' 2 {x,y) mSE ) from f c {po) mSE , we find the 2-point distribution 

t 2 (po) mSE = r c (p r SB -r 2 (po) mSi3 

4 /*°°, 9(s-m 2 ) fl 3m 2 , m 4 l . A/ s (4.11) 



As a next task, we compute the integral of the principal value part of ( 4.11] ), 
denoted by X(po): 

w , -ir 4 _ p _ 2s 2 - 3m 2 s + m 4 
X { p ) = —p Pj m ds s3{s _ p ^ 



4tt 7 m 2 [ s s 2 s — Pq 
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with 



a :- 



3m 
Po 



fi 2 ' P ' 

Po Po 



3m 
Po 



m 
Vo 



7 : = 



m 
Po 



Integration of the partial fractions in ( 4.12j ) yields 

,2 _4 



X(p ) 



iT 
2^ 



2 3m" m 
P °-^- + 2^' l0g 



2 2 

pg — m 



rm 



+ 



m 



5p| 
4 



(4.13) 



(4.14) 



4.2 Matter Self-Energy Two-Point Distribution and Freedom in 
its Normalization 

From Q) with (g3| ) and (PD , we can derive the 2-point distribution for the 
matter self-energy 



inSE 



(4.15) 



and in an arbitrary Lorentz system, the matter self-energy distribution for p £ 
V+ is 



r 

2^ 



3m 



m 

"2~ + V 



log 



2 2 

p — m 



z7T0(p 2 — m 2 ) 



+ 



m 
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5p 
(4.16) 



The obtained loop contribution is UV finite and cutoff-free due to the causal 
scheme. 

This loop contribution was also calculated in [22] within the operator regular- 
ization scheme. Parts of their result agree with our expression in ( 4.16Q , whereas 
differences concern the explicit presence of parameters, of other p-dependent log- 
arithms and terms which go as p in their expression for S(p). In the causal 
scheme, these latter cannot appear, because the singular order remains the same 
after distribution splitting. 

The retarded part in (^^) is the boundary value of the analytic function of 
complex momentum p + iff, r\ = (e, 0), e > 0: 



f(p) c 



iT 
2^ 



9 3m m 



2p 2 < 1 ° g < 1 



2 \ 2 

m \ m 



5p 2 



(4.17) 



Having split d™ SE with uo = 2, an ambiguity in the normalization of E(p) of the 
type 



N(p) = — (c + c 2 p 2 ) 



(4.18) 



must be taken into account. In order to fix the constants cq and C2, radiative 
corrections to the matter Feynman propagator by matter self-energy loops are 
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considered: letting formally g — ► 1, they are of the form 



i D^(x - y) + jd 4 xid A X2 ( - i D^{x ~ x i)) 



x i [E(xj - x 2 ) + iV(xi - s 2 )] ( - i D F Jx2 ~ y)) + ■ ■ ■ ■ (4.19) 
In momentum space the series becomes 

D*Sp) + DM (2^) 4 (E(p) + Nip)) D F m (p) + b F m (p) x (4.20) 
x (2^) 4 [±{p) + Nip)) DUP) (2^) 4 {t{p) + Nip)) D F m (p) + ...=: E(p) tot , 



where D^ip) = (2ir) 2 (—p 2 + to 2 — iO) 1 . The geometric series in ( |4.20D leads 
to 



E(p) tot = DUP) (1 + (2^) 4 (S(p) + 2V(p)) E(p) 
1 1 



tot 



(4.21) 



(2vr) 2 - p 2 + m 2_ i0 _ (2vr)2 (fj(p) + JV(p)) 



Mass corrections are avoided by requiring 

[t(p) + Nip)] 

This implies 



2 / 3 
(•(, m ( - - c 2 



(4.22) 



(4.23) 



so that the matter self-energy distribution including its normalization reads 



Zip) + N(p) = — 



3m 2 



2 + 2p 2 ' 



log 



2 2 

p — m 



77 V - 



i 7T 0(p 2 — m 2 ) 



+ ip 2 - m 2 ) (pa - -) 



+ 
(4.24) 



For the remaining constant c 2 , we cannot provide here a value. A restriction 
should come from the vertex corrections A{p, q) a p to third order in the 3-point 
distribution Ts(x, y, z) =: 0(x)0(y)/i a/3 (z) : h{x,y,z) a p. 



5 Perturbative Gauge Invariance to Second Order in 
the Tree Graph Sector 

In this section we show that the condition of perturbative gauge invariance to 
second order in the tree graph sector generates a local quartic interaction of the 
form ~ k 2 : hhcfxp : 5 (4) (x — y), which agrees with the second order C M in the 
expansion of the matter Lagrangian density C M , Eq. fl2.24|) . 



26 



5.1 Methodology 

Since R 2 (x, y) is trivially gauge invariant due to its definition and to the gauge 
invariance of Ti(x), instead of Eq. (2.20) for n = 2, we have to examine whether 

d Q R 2 {x } y) + d Q N 2 {x,y) = dlR v 2/l {x,y) + d y v R v 2/2 {x,y) + 

+ W /1 (x,y)+W /2 (s,y) (5.1) 

can be satisfied by a suitable choice of the free constants in the normalization 
terms N 2 , N 2 / t and N%, 2 of the retarded parts R 2 , R 2 h and R 2 / 2 - Here, R 2 ^ 
and i?2/2 are ^ ne re tarded distributions obtained by splitting the inductively 
constructed distributions 

D» 2/1 (x,y) := [T^ixlT^y)] and D» 2/2 (x,y) := [T ± (x) , (y)] . (5.2) 



This procedure has been described in |10[] for pure QG. It turned out that 
Eq. (5A) can be spoiled by terms with point support ~ <5 (4) (x — y), only. N 2 , 
N^h and N%/ 2 are by definition local terms, but there is another source of 'local 
anomalies', namely the splitting procedure for tree graphs: in the inductive 
construction of R 2 h, the Q- vertex Ty x can give rise to expressions of the form 

D v 2/1 {x,y) =:0{x,y): d u x D m (x - y) . (5.3) 

:0(x, y) : is a normally ordered product of four fields, being the other two fields 
in the commutator (|5.2;). m can also be zero, if the commutator contains two 
graviton fields. The retarded part is simply 

I% n (x,y) =:0(x,y): %D%(x-y), (5.4) 

because the trivial splitting for tree graphs follows from D m = — with 
the correct support properties as explained in Sec. |2.1| . Applying the derivative 
dy, which forms the divergence in d^R^n, we get 

:0(x,y): a x D^(x-y)=:0(x,y): [5^ (x - y) - m 2 D%] . (5.5) 

The expression : 0(x, y) : <5 (4) (x — y) is a local anomaly. A corresponding mech- 
anism works also for R 2 / 2 - 

We define by an(<9jj;i?2/i + ^vR 2 i 2 ) the set of all local anomalies generated 
by the described mechanism. 

Following [10 1, gauge invariance is preserved if we can choose N 2 , N 2 / 1 and 



^2/2 so that the condition 



d Q N 2 (x, y) = Md X uR V 2/i + d y v R u 2/2 ) + W A (x, y) + W /2 (x, y) (5.6) 
involving the local terms of ( |5.1|) is satisfied. 
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Note that in QG coupled to matter cLqR2 does not generate local terms with 
matter fields involved, namely of the type : uhcjxj) :. This is in contrast to the 
much more involved pure QG case. 

At this point we realize that it is not sufficient to consider Tf 1 only. Also the 
graviton and ghost first order couplings and T" have to be taken into account, 
because they yield also local anomalies with external operators ~:uh(j)<j>: when 
splitting the commutators of ( |5.2| ). 

Using a simplified notation which keeps track of the structure of the coupling 
only, gauge invariance to first order then becomes 

d Q (:hhh: + :uhu: + :h<f><t>:) = d£ ( : {uhh}" : + :{uuu} v : + : {u(/)(f)} v :) . (5.7) 

s v / s v ' 

rph+u+M mf h+u , Tiy M 

With the 'extended' Q- vertex Ty^ +U + T",™, we decompose the commutator 
defining -Dg/i m ^ ne following way: 

Z?2/i(£,y) = pure QG sector + [ :{uhh} u :,:hcj)cf): ] + 
+ [:{«#}" :,:«/»«:] + [:{t ' "' ' 



In the pure QG sector, which involves terms of the type : uhhh : and : uuuh :, 



perturbative gauge invariance has been shown in [10|. 



5.2 Explicit Calculations 

Our task consists in the investigation of the three remaining sectors in Eq. ( |5.S| ) 
in which matter and graviton fields are mixed together. We denote the three 
commutators in ( |5.S| ) as the graviton-, ghost-, and matter sector, respectively. 
Performing one contraction they lead to 

Dtf 1 (x,y) = (:uh,^: [h(x) , %)] ) " , 

X y 

D Z/l(x,y) = ( -.^h^- [u(x),u(y)]) u , , 5 9 ^ 
x v 

D^(x,y) = {: u4> hep : [0(x), 0(y)] ) v , 



respectively. By considering the explicit form of Tf in [1C] and Ty™ in ( 2.19| ), 
we find that no local anomalies arise in the ghost sector. 

Let us compute the local anomalies in the graviton sector. From the ex- 
pression for Ty 1 +U in [10] we isolate the terms that generate local anomalies 
according to the mechanism described in Sec. [O]. They are 

T-h+u = f + 1 :u i* h P* h /»: + l -u^h?*: -- :u^hhy. 

1 i , (5-10) 

- - :u»h tll h tV : +- :u%h» p h v : - \u%bP°h^ :+...}. 
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Then, D%k contains the following terms that generate anomalies 

2 2 

:u"/i iM #: -2 :v? p hr$y$y. +m 2 :/ p /i w #: } fl£.Do(a: - y) • (5.11) 



m 
IT 



The first two fields in the normally ordered products depend on x, whereas the 
matter fields depend on y. The retarded part R^/i nas exac tly the same form 
as in ( 5. 11| ) with the replacement Dq — > -Dg et . Applying <9J to R%j[ (5.11), we 
obtain the local anomalies (local anomalies coming from d^R^k are just the 
same, therefore we get a factor two): 



an(3*i?^+d^; 2 ) = y :,/';,//'"- o,.o.,: + :<'" //'>,.. o,, : + 



m 

T 



??? 



— -2 lu^hP^^: +m 2 :uf p /^#: j - y) . (5.12) 

Because of the <5-function, all fields have now the same space-time dependence. 

In the matter sector, the first term in the matter Q-vertex of Eq. ( [2.19D is 
the only one that can generate local anomalies, because it carries the z/-index 
as a derivative. Then, Dyf becomes 

D^{x,y) = -^-{:u^^cj>y d* + ^ :u»<t>yh<j>: } S£ D m (x - y) . 

(5.13) 

The contribution coming from has x «-> y and the derivative attached to 

the first term is d„ . The retarded part R^Jf nas the same form as in (fyJD with 
the replacement Dq — ► -Dq 6 *- Since 

:A(x)S(!/): ^<*>(x - y)+ : A(y)5(x) : dl^\x-y) = 

= - :A(x) >a B(x): 5 w (x - y)+ : A(x)B(x) j0l : 5 (4) {x-y), (5.14) 

by applying d% to R v 2 f[ and dy to R%j%, we obtain 



2 r 



MdZRyl + W/2 ) = y [ - ™ 2 -^Wy4>- + + 

+ :u tk hP a 4>^' ~ ■u^h^^y. - :u^6 pa (t>y] 5 (i \x - y) . (5.15) 



Because of the 5-function, all fields have now the same space-time dependence 
and have been recast in the form :uf 
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5.3 Quartic Normalization Terms 

Summing up, we have found all the local anomalies arising from d^R^H+dvR^ 
attached to normally ordered products of the type : uhcjxft :. These can be 
organized in three different types according to their Lorentz structure apart 
from derivatives: type I :u fJj h f " T 4>(p: with derivatives dp,d p ,d a ; type II :u^h(j)(j): 
with derivatives dp,d p ,d p and type III : u^W^ : with derivatives d p ,d a ,d a . 
Different Lorentz types do not interfere. 

From Eqs. ( |5.12| ) and ( |5.15 ), the local anomalies of type I are 



an(^^ A + d y v R. 



v ■> 
2/2/ 



type I 



K 
Y 



: u 



■■ + 



5 {A) {x-y). (5.16) 



From Eqs. ( |5.12| ) and ( p. 15 ), the local anomalies of type II are 
an(^ /1+ ^ /2 ' 



type II 



and those of type III are simply 



K 

y 



2 ^2 

:u^h<j)(f): +— :u^h 



2 '>>* ~ T T ' ' 2 
+ m 2 -.u^hcj) u4>'- 



> : + 

5 w (x-y), (5.17) 



an(^/i + W /2 



/2 



type III 



m 2 :vr p h w 44\ 5 w (x-y). (5.18) 



According to Eq. ( p.6[ ), the question is whether these local anomalies can 
be written as divergences and therefore can be compensated by corresponding 
local divergence terms coming from d^Ny^ + d„Ny 2 . If it is n °t possible to 
reach such a compensation, normalization terms on the left side of ( |5.6| ) have to 
be introduced in order to preserve gauge invariance. 

Due to the identity 

d x y (:A{x)B{y): 5 w (x — y)) + <9^( :A(y)B(x) : 5 w (x-y)) = 

=:A(x), u B(x) : 5 (4) (x - y)+ :A(x)B(x) >v : 5 (4) {x - y) , (5.19) 



the local anomalies of type I can be indeed written as a divergence: 
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The same is true for the local anomalies of type II: (5-21) 
Eq. ([5T7D -^m 2 d*(:u^h 4>4> : 5 w (x-y)) +m 2 d£( : u^h # : 5 (4) (x-y)) . 



^ y y x 

Therefore, by choosing the appropriate local divergence normalization terms 



d^N^ft + d„Ny 2 , Eq. (5.6) restricted to the Lorentz type I and II holds. 



The local anomaly of type III, Eq. ( 5.18| ) cannot be written as a divergence. 



Eq. ( |5.6| ) forces us to consider normalization terms N2- Then, we require that 
the local anomaly of type III has to be coboundary, namely a gauge variation 
of a normalization term N2: 

2 2 

d Q N 2 (x,y) = :u^{x)^{x)4>{x)4>(x): 5< 4 \x - y) . (5.22) 

The only possible N 2 that satisfies this requirement is 

N 2 (x,y) = -\ :hh^: } 6™(x - y) . (5.23) 

Taking the factor \ for the second order iS-matrix expansion ( |2.1| ) into account, 
the quartic interaction in ( |5.23| ), quadratic in k, generated by gauge invariance 
agrees exactly with the term of order k 2 in the expansion of the matter La- 
grangian density L (HH). But in our case, this mechanism of generation 



works in a purely quantum framework. 

The only objection to this result is that this N2 is not a 'proper' normaliza- 
tion of a tree graph in T 2 , obtained starting with T^ +u and T™. But it can be 
considered as a normalization term of box-graphs in fourth order with external 
legs ~:hh(f)(f): constructed with two T^ +u and two T X A/ or with four . 

5.4 Massless Matter Case 

For the massless matter coupling T x M (x) = i| :h a ^(p ia ^)^:, perturbative gauge 
invariance to first order reads 

d Q Tf(x) = dl(^ { -.uP^y. -~ :u», p : }) =: d x u T^{x) , (5.24) 

namely, the matter Q-vertex of ( |2.19| ) in the m — ► limit. Performing the same 
calculation as for the massive case, it turns out that the local anomalies are those 
of type I only, Eq. ([5 . 1 6|) which can be written as a divergence and therefore 



compensated by proper normalization terms d^N^^ + d^N^^ as i n ( 5.20 ). From 
Eq. (|5.6|), it follows that dqN 2 = 0, which is certainly satisfied by N2 = 0. This 
agrees with the fact that even classically there are no hhcfxfi couplings in the 
m = case: the expansion of C M written in terms of the Goldberg variable 
reads 

£m = \<t>, P <l>' p + 1 h^M* ■ (5-25) 

Therefore, Cff = 0. This concludes our discussion of the condition of perturba- 
tive gauge invariance to second order for tree graphs. 
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Appendix 1: The I^(p)... -Integrals 

The I^{p)... -integrals are defined by 

^W-Mr/w/rP* ■= *«P ~ «)* -m 2 )0(± (p° - q )) 5(q 2 - m 2 ) 

x @(±q°) [1, q P , q»qv, q^q v q P , q^,q v q P qa\ ■ 

(A.l) 

Let us investigate 1^ (p) in detail, because it contains the main feature of the 
calculation of the other integrals, p is time-like due to the presence of the 0- 
and ^-distributions in the integrand. Therefore, we may choose a Lorentz frame 
such that p v = (j>o, 0), po > 0, then 

4 +> (n>) = /A«(rf - 2p * + i, - l</l 2 - ™ 2 ) &(p« - m)S(<! 2 - ™ 2 ) e(+») 



3 5(p 2 - 2p E q ) G(p - E q ) 



2E q 



(A.2) 



with E q = qo = \J q 1 + m 2 . Then 

/;:'(Po) = 4,pi,i^i a (|- E ,)e (M -£ ( 

PO Jo E q \ 2 



(A.3) 



because of the (^-distribution 



- v 7 ^ + /m 2 = ^=H?| = ;/^-m 2 . (A.4) 



Therefore 



4 +) (Po) = - B(p 2 - 4m 2 ) 6(p ) M-m 2 
P0 V 4 



7T 



- -^/l - 1^! 6(p 2 - 4m 2 ) 9(p ) 
Po 



(A.5) 
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E m ' can be calculated analogously and the result in an arbitrary Lorentz frame 
reads 



I™(p) = \ Jl - ^ 9(p 2 - 4m 2 ) 0(±po) =: I d(p)W . (A.6) 

Computing I^ip)^ for = (pojO), po > 0) we have a non- vanishing con- 
tribution only for fj, = 0. An additional factor go is therefore present in the 
integrand in ( |A.3j ) , which is set equal to E q and then to ^ . This leads to 




IL ±} (P), = P» t \i l ~ ^ " W ) (±Po) = £p„ d(p)£ } , (A.7) 

in an arbitrary Lorentz frame. 

For I^(p)[_iv, we have to take the two conditions 

f 4 ±) (p)/ = "» a 4 ±) (p). fA8) 

into account. The first condition is imposed by the presence of 5(q 2 — m 2 ) 
in ( |A.1| ), Computing p^p u I^(p)p, v for Pu = (?>o> 0), po > 0, we get an additional 
factor (poE q ) 2 in the integrand ( [A .3] ) which gives ^ due to — E q ). In a 
general Lorentz frame, the second condition in QA.8D must hold. Therefore, 
making the ansatz 

#° (pW = Hp 2 ) p»p» + &(p 2 ) p 2 Ir) / i ±) (p) , (A.9) 

the conditions ( |A.8[) imply 

«« = 5-$, *» = i* + £. (A.10) 



so that 



(A.12) 



For I^{p)aupi we have to take the two conditions 
f I^(p), v v = m 2 I^(p),, 

into account. The first condition follows from the definitions of the /^-integrals. 
Computing p^p u p p I ( r ^ ) (p)^ p for p v = (po,0), Po > 0, we obtain an additional 
factor (poE q ) 3 in the integrand ( |A.3 ) which yields ^p. Making the ansatz 

4 ±} 0W = (c{p 2 )p^p u Pp + d{p 2 )p 2 (p p r]nu + PpVpv + PuVph)) I^Kp) i 

(A.13) 
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then ( |A.12| ) imply 

4 2p 2 ' yy ' 24 6p 2 

and therefore 



<P ) = 7 ~ 753 > c (p ) = 777 + 7^2 > ( A - 14 ) 



7T 



(A.15) 

For I { ^ ) {p)p LV pa- 1 we repeat this calculational scheme again. For the reasons 
pointed out above, now three conditions 

' I^(pW P = m a jnW(p) /lw> 

< I^(pW = f^I^(p), (A.16) 

k pVpV / w (p W = g (p) , 

must hold. Therefore, the ansatz 



I m\p)^P^ = (e{p 2 ) PpPuPpPa + f{p 2 )p 2 {PpPaT]pu + PpPpVau + PpPuV^fi + 

+PpPaV P u+PuPar]pp+PpPur]pa) + g(p 2 ) P 4 " {r]p p r]ua + VpaVup + VpuVpa)) ^i ±] (p) , 



(A.17) 



leads to 



, Q N 1 3m 2 m 4 „. 9v —1 7m 2 m 4 

<P ) = 7 " TZ^ + O ' /(P ) = IK + 



5 5p 2 5p 4 ' ^ VjK ' 40 60p 2 15p 4 
1 / 1 m 2 
15 V 16 ~ 2p2 + 



, 1 / 1 m 2 m 4 i 

9(P 2 ) = 77 T7.-^ + -T > (A-18) 



so that 



4 '(pW = 1 1 1 - ~pT + ZJ ) P»P»P P P° " ^ " TV + V J P 

X {PpPaVpu +PpPp,Vav +PpPvVvp +P^PaVpu + PuPaVpp + PpPuVpa) 

(A.19) 
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Appendix 2: The J (±) (p)„. -Integrals 

The products of Jordan-Pauli distributions of Eq. (4.2) are evaluated in mo- 
mentum space: 

<5f (p) = J^y> Jd\b^(p- q ) ft* (q) 

4 jd 4 q5((p- q ) 2 )Q{±(p - qo ))5(q 2 -m 2 )G(± qo ) (B.l) 



-1 

: W) 1 



and 



^(P) = t^ti ld\b^{p -q){-iq a ) D%{q) 



+i 

(27T) 4 
+* 



(27T) 



y*d 4 g g a J((p - g) 2 ) 9( ± (po - go)) % 2 - ™ 2 ) e(±go) (B.2) 
4 / (±) (p)*- 



Let us calculate / (+) (p). From the 5- and ©-distributions, it follows that p is 
time-like. We choose a Lorentz frame in which p u = (po, 0), po > 0. Then 



! {+) (po) = Jd A q5{pl - 2p q + q 2 - q 2 ) B(p - q ) e(g ) % ° ^ 

= jSr 6{p2 ° ~ 2poEq + m2) @{po ~ Eq) ■ 



2Eq <-> 



With E q = \J q 2 + m 2 . The ^-distribution implies po = E q + \q\. From E 2 = 
q 2 + m 2 = (po — \q\) 2 , we obtain \q\ = (p 2 , — m 2 )/2po, which, due to po > 0, 
yields 0(pg — w 2 ) and 

r°o |„|2 

/( + )(po) = 2vre(p 2 -m 2 )e(p ) / 5(p 2 - 2p 0J E q + m 2 ) 

= 2vre(p 2 -m 2 )e(p ) /°dE, M S(E q 



9 9 

Po + m , 



2p 



^e(p 2 -m 2 )e(p ) A /^- 



E = a, +m 



(B.4) 



= * e(p 2 - m 2 ) e(p ) (p§ - m 2 ) = | e(p 2 - m 2 ) e(po) (i - ^) . 

zp z p 

Therefore, in an arbitrary Lorentz frame: 

/W(p) = | (p 2 - m 2 ) G(±p ) (l - ^) • ( B - 5 ) 
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The second integral can be in a similar manner computed: in the Lorentz frame 
with p u = (po,Q), po > 0, I w (po)i, i = 1,2,3, vanishes for symmetry reasons, 
then 



I w (po)o = E q 5{pl - 2p E q + m 2 ) 6(p - E q ) 

roo I |2 

= 2^e(^-m 2 )6(p ) / d\q\ 7 P- r E q 6(E q - 



pfi + m z ^ 



2po 

= JL e{p 2_ m 2 )e{po) r dEqlqlEq6{Eq _Po+^ ) (B . 6) 

PO Jra ^P0 

= _ efa ,- m ) e( „ ) )_ 5 _ 

= ^e«-m 2 )6 W (l-^) (l + ^)- 
Therefore, in an arbitrary Lorentz frame we have 

I {± \p) a = \pa U ~ ^) 0(p 2 - m 2 ) 9(±po) . (B.7) 
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